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Flavor symmetries successfully explain lepton and quark masses and mixings yet it is 
usually hard to distinguish different models that predict the same mixing angles. Further 
experimental input could be available, if the agents of flavor breaking are sufficiently low in 
mass and detectable or if new physics with non-trivial flavor charges is sufficiently low in 
mass and detectable. The recent hint for lepton-nonuniversality in the ratio of branching 
fractions B —>■ I\ pp over B —> Kee, Rk, suggests the latter, at least for indirect detection 
via rare decays. We demonstrate the discriminating power of the rare decay data on flavor 
model building taking into account viable leptonic mixings and show how correlations with 
other observables exist in leptoquark models. We give expectations for branching ratios 
B —> Kit', £?( s ) —> It! and t —> t! 7, and Higgs decays h — > IP. 


I. INTRODUCTION 

The LHCb collaboration [T] has presented recent data on lepton-nonnniversality (LNU): Rk, 
the ratio of branching fractions B —> ; over B —> Kee, Rk has been measured as 

R\ flCb = 0.745 ±g;{$g ±0.036 (1) 

in the dilepton invariant mass squared bin 1 GeV 2 < q 2 < 6GeV 2 , which represents a deviation of 
2.6(7 from lepton universality [2J. In view of its generic leptoquark interpretation [3] we investigate 
directions in flavor model building that link the electron-muon non-universality to other flavor 
observables particularly lepton masses and mixing and lepton-flavor violation (LFV). In general if 
LNU is present, LFV can be expected, too [4J. Rk has also been addressed within (composite) 
leptoquark models j5]; they do not aim at explaining LFV. Also recently, [6] considered scalar 
leptoquark contributions to rare leptonic B meson decays. 
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Here we consider frameworks where both phenomena (LNU, LFV) stem from the same (scalar) 
sector. To be more specific, we investigate the possibility of simultaneously accommodating LNU 
as observed through Rk, and LFV, as observed in neutrino mixing and as hinted at 2.5 a by CMS 
in h—^Tfi (7]. We do so by adapting frameworks with non-Abelian flavor symmetries that predict 
the leptonic mixing matrix. 

Due to the Standard Model (SM) gauge group, the maximal flavor symmetry is an 17(3) for each 
chiral quark and lepton species. This observation motivates the use of SU{3)f or subgroups. After 
discussing briefly a framework of leptoquarks within the continuous group SU(3)f, [E] we consider 
explicit realizations with the discrete group A 4 in frameworks within supersymmetry (SUSY) [5fU2|. 
While the specific results are related to the details of the models, our general conclusion is generic 
even if other discrete subgroups of SU(3)f are invoked. If the same flavor symmetry breaking fields, 
often referred to as flavons, are involved in Higgs and leptoquark couplings, the respective Yukawa 
structures are related and there are links between what is predicted for LNU and LFV (in lepton 
mixing and Higgs decays). 

The paper is organized as follows: in Section [TT] we identify viable, data-driven patterns for flavor 
structure of the leptoquark to fermions coupling matrix. Predictions for rare LFV decays of leptons 


and H-mesons are given in Section III The results of these sections are rather independent of flavor 


models. In Section IV these benchmarks are compared to what is realized in theoretical flavor 
models, where we consider models based on SU(3)f and A 4 known for explaining lepton flavor. 
Rare Higgs decays are addressed in Section [V] In Section VI we conclude. In several appendices we 
give formulae and auxiliary information. 


II. DATA-DRIVEN FLAVOR PATTERNS 

We consider models with a scalar 51/(2) doublet leptoquark A(3,2) 1/ / 6 , with electric charges 
(2/3,—1/3) and one with a scalar SU( 2) triplet leptoquark A(3,3)_ 1 / 3 , with electric charges 
(2/3, —1/3, —4/3). Such frameworks correspond to the RL model and the LL model, respectively, 
of [3], which discussed them as “single particle” explanations of Rk • For this reason, apart from 
having renamed the leptoquarks as A, throughout this section we use the notation from [3]. 

In view of the uncertainties in the current data in the following we neglect factor of 1 /2 differences 
in the semileptonic 4-fermion operators after Fierz-ordering and other small differences between 
the RL and LL model variants. This way the phenomenological constraints on the leptoquark 
parameters relevant to our study give identical hierarchies in flavor patterns, which we hence study 
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jointly. We stress that despite the present day patterns being similar, there are differences in terms 
of model-building. The 517(2)-doublet leptoquark model couples to right-handed quarks 1 

£ = -A q t A ( qPL£) + h.c. , (2) 

while the SU (2)-triplet leptoquark model couples to left-handed quarks 


C = -\ qe A* ( q*i) + h.c . . 


( 3 ) 


In each case, the couplings to the leptons are left-chiral and quark (q) and lepton (£) denote flavor 
(generation) indices. We denote the mass of the leptoquark by M. 

We continue discussing constraints on the generalized Yukawa A q i for q = d,s,b and £ = e,/ a, t: 


Rk implies at 1 a [3f: 


A = 


A de A dfj, A dr 

Ase A sm X st 
y A be Aft^j Xb T J 


0.7 <Re[A se Aj* e - X sll X* b / U ™ )2 < 1.5 , 


( 4 ) 


( 5 ) 


We restrict our analysis to M > 1 TeV to be conservative with collider bounds m-nn, hence, in 
view of Eq. 0, |A se A^ e - X sll X* bfl \ > 2 • 10 -3 . 

There is also a contribution to B s mixing from a box diagram with leptons and leptoquarks in the 
loop giving rise to an operator of the form bsbs with the complex coefficient A S ^A^) 2 /(167r 2 Af 2 ). 
Numerically, 


(AseAfte + X s ,X* bfl + X ST X l T Y 


M 2 


(12 TeV) 2 ’ \ (!7.3 TeV) 2 


M 2 


( 6 ) 


where the weaker bound is obtained from the mass difference |Am^ p /Ar?Tg M | < 0.15, whereas the 
stronger one in parentheses stems from the upper bound on the B s — B s mixing phase —0.015T0.035 
(defined relative to the SM phase in the b —> ccs decay amplitude) [18]. The complementarity of 
the | Ail| = 1 and |A1?| = 2 constraints allows to fix the mass scale and the dimensionless couplings 
separately, yielding an upper limit M < 0(50) TeV. Couplings to different quark flavors are within 
four orders of magnitude, 10 4 > X s i/Xu > 10 -4 [3]. 


1 The RL model has similarities with the R.-parity violating minimal supersymmetric standard model (MSSM); 
the corresponding term is the A 'LQd c coupling. Note also that the A(3, 2)\/e with mass much lighter than the 
GUT-scale does not require further model-building to avoid proton decay m ■ 
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The simplest scenarios are when the leptoquark couples to one generation of leptons only, as 
considered in (3]. In this case A reads, 

/ Aj, D n \ ^0 (A 

AM = 




A de 0 0 
A se 0 0 
\ Afee 0 0 J 


0 \ S fj, 0 

0 A bp 0 j 


( 7 ) 


A third variant with couplings to r only cannot explain non-universality between electrons and 
muons as hinted currently by Rk data. With either of those limits, A^l or A^, there is no LFV in 
rare processes. 2 

Generalizing the single lepton scenario we consider the hierarchical pattern 


AW ~ A 0 


( Pdn Pd Pd ^ 


PK p p 
K 1 1 


( 8 ) 


Here, we allow for quark flavor suppressions pd = Xdi/Xbl and p = X s £/Ay having larger couplings 
for heavier generations. This is suggested by the observed quark mass pattern; see Appendix [B] 
for details on how this is realized with the Froggatt-Nielsen (FN) mechanism |20| from [/(l)-flavor 
symmetries. Here we simply take typical phenomenological profiles for Higgs-Yukawas of quark 
doublets (q) and down-type singlets (d) from [2XJ, and apply them to the leptoquark-Yukawa, 
leading to 


p ~ e 2 , p d ~ e 3 or e 4 , (q) p ~ e, p d ~ e or e 2 , ( d ), 


(9) 


where e denotes a flavor parameter of the size of the sine of the Cabibbo angle, e ~ 0.2. We stress 
that what we mean by the hierarchies and flavor patterns in general is that they hold up to numbers 
(including CP phases) of order one. The overall scale Ao in our ansatz (A^) depends on M and is 
fixed by Rk through Eq. 3 


2 Generically Rk < 1 can be explained by a suppressed branching ratio to muons, or an enhanced one to electrons, or 
combinations thereof, the details of which driven by an underlying model of flavor. Present data, however, is much 
more copious for rare decays into muons than into electrons due to the large data samples from the hadron machines 
at Fermilab and CERN LHC. This is likely to change in the nearer term future with the Belle II experiment taking 
data. As a consequence at present the constraints on the b —> see modes are much more loose. In addition, hints 
for new physics in B —> K*pp angular distributions could point to a simultaneous explanation with Rk, see m 
and references therein. So, the rather minimal scenario where Rk originates at least predominantly from muon 
mode suppression and the electrons are SM-like presently has phenomenological support. 

3 Note that here we do not explicitly specify U(1)fn charges for the leptons i, whose assignments are not fixed 

by phenomenology as they depend on the model of neutrino masses. While in general the texture Eq.|ij| with 
suppressions Eq.([9j is expected to hold, there is the possibility for cancellations between quark and lepton charges, 
that could lead to inverted hierarchies in the leptoquark Yukawa A. We do not consider such solutions; they are 
absent as long as all charges are non-negative. 
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observable 

current 90 % CL limit 

constraint 

H(/x -» ei) 

B(t —> eq) 

B(t -» iij) 

B(t ->• 

B(B -> K^e?) 
B(B -4 Jir ± e T ) 
8(B->A> ± r :F ) 
B(B — > 7r/.< ± e =F ) 

5.7 ■ IO" 13 P2 

1.2 • IO" 7 PH a 

4.4- 10" 8 p5] 

6.5- 10" 8 P3 

3.8- io- 8 pg 

3.0 ■ IO" 5 PI 

4.8 • IO" 5 |H] 

9.2 • 10" 8 PH 

1 \ \* 1 M 2 

kV'V/L ~ (34TeV) 2 

1 \ \* 1 M 2 

l / '?e / VI ~ (0.6TeV) 2 

1 \ \* 1 ^ M 2 

\ A q^ A qr\ ~ (0.7TeV) 2 

\\ \* | <- M 2 

1 a s[i a st 1 ~ (3.7 TeV) 2 

a/ A SjU A be + A& M A se ^ (19.4 TeV) 2 

\/|As T A 6e + Ab r A se < (3.3TeV) 2 

y/l As,u. A 6i _ | + |Ab M A ST | (2.9 TeV) 2 

x/IAd/xAjJ + A6 M A de ^ ( 15 . 6TeV ) 2 


future sensitivity 

6 • io -14 pg 

5 • io -9 m 
2- io- 9 pg 


“In previous versions of this work this bound was incorrectly stated as 1.2 • 10 8 , which reflected the abstract of 
the original version of |24j in the arXiv, which has been corrected in June 2015. 


TABLE I: Selected LFV data, constraints and future sensitivities. Here, q = d,s, b. The Belle II projections 
*2(1! are for 50 a& _1 . For the constraint from B(t —» /i?y) we ignored the possibility of cancellations with 
see e.g., |30'. Following [3T] we ignore tuning between leading order diagrams in the £ —> £ >r y 
amplitudes. 


We emphasize that our main goal is to simultaneously explain lepton mixing and LNU. Lepton 
mixing is structurally very different from the hierarchical quark mixing matrix, which suggests a 
different origin of flavor for the leptons. We therefore consider in Section |TV| scenarios where the 
entries between columns of A are related by factors that are ±1 or 0 as in Eq. ([7]). Such structures 
can occur naturally from a non-Abelian flavor symmetry distinguishing the generations of leptons. 
We state the set of assumptions that lead to an effective factorization of quark flavor (rows) and 
lepton flavor (columns) in A: the Higgs, leptoquark and quarks are singlets under the leptonic flavor 
symmetries while leptons are neutral under the quark symmetries. Most scenarios we consider will 
be factorized, including the hierarchical ansatz Eq. ([8|. Therein the flavor parameter k allows 
to split the electrons from the other leptons, and is further discussed together with LFV below. 
Although it is not completely general, the structure of A^ is useful to obtain predictions, which 


we do in Section III that can apply for more general flavor patterns, as shown in Section IV D 
Non-factorized scenarios are discussed in Section |IV C 1| and Appendices [D] [E] 

If couplings to more than one lepton flavor are present, existing constraints on LFV processes 
need to be considered. In Table [I] we collect LFV data, the corresponding constraints on A-matrix 
elements, and give future sensitivities whenever available. Constraints on couplings to r-leptons 
exist from the H-factory experiments Belle and Babar [211I25J [27]. These bounds are not competitive 
with 6-physics ones but can be improved in the near term future at Belle II |26| . The strongest 
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constraint on X qT is therefore from £> s -mixing Eq. Q. Consequently, phenomenology does not 
currently require a split (within our approximations) between muons and taus, which is why we 
refrain from adding a flavor factor separating the second and the third column of A^ in Eq. |8]). 
Note also that while it would be straightforward to implement such a factor, for leptoquark masses 
near maximal, Rk already requires muon couplings X qfl of order one, cf. Eq. (J5|, so allowing A 
significantly larger than X qfl would challenge perturbativity. 

The strongest constraint on LFV stems from p, —> eq from MEG 2013 

M 2 


qr 


IVA^| < 


(34TeV) 2 


q = d,s,b. 


( 10 ) 


This is only mildly stronger than the Rk bound Eq. ([5]) , implying k/ p < 0.5. The requirement 
of flavor hierarchies between electrons and muons can become more pronounced in the short-term 
future, with MEG upgrading their sensitivity to 6 • 10 -14 [25] , which implies (k/ p) U pgrade 1$ 0-2, 
if no signal is seen and Rk remains close to its current determination. If the // —>• e 7 constraint 
tightens, it implies that in Rk there has to be a dominant contribution to one lepton flavor with 
one branching ratio closer to its respective SM branching ratio than the other - ruling out solutions 
to Eq. ([ 5 ]) where e and p couplings are of roughly the same order of magnitude. 

Current limits on LFV in H-decays B(B —> KR^l'^) are also given in Table [j] They are weaker 
than Eq. ([ 5 ]). Comparison of the bound on B —> Kp^e^ with the one on B —> np^e ^, also shown, 
implies p d /p< 1.6. 

Rare kaon decay data provide the strongest constraint on couplings to d-quarks [31] : 


| A dfi A 


I rS 


M 2 


(183TeV) 2 ' 


( 11 ) 


Due to chiral symmetry contributions from purely left-chiral leptons to the anomalous magnetic 
moment of the electron or the muons A~ |A 9 z| 2 /(167r 2 )m 2 /A/ 2 are much smaller than current 
experimental sensitivities. 

To summarize, the phenomenologically viable range for A^l parameters in Eq. ([8]) is 


Pd <0.02, k< 0.5, 10” 4 < p < 1 , n/p< 0.5, p d /p< 1.6. 


( 12 ) 


In view of p —> eq, if p <C 1 (corresponding to a low mass leptoquark), then necessarily k< 1. 
Regarding p —> e conversion in nuclei, the current best limit T(pAu —> eAu)/T(pAu) < 7 ■ 10 -13 
is consistent with the bounds (12), however, future experiments such as COMET [33| and Mu2e 
with sensitivity below ICE 16 are sensitive to the parameter space of leptoquark models discussed 
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here 4 . We note that corresponding tree-level contributions can always be avoided by decoupling 
down-quarks, suppressing pp, or decoupling electrons, suppressing k. 

Note that the model by [5] with the leptoquark being an SU(2)l triplet has a flavor structure 
Xij ~ e n(?d e ^(«) ; where the e-powers are adjusted to give the correct fermion masses, similar to what 
could be obtained through a FN mechanism with [/(l)-flavor symmetries. The resulting pattern 
turns out to be a viable subset of the one presented in Eq. ([8]) . Specifically, p and pd follow Eq. (|9| , 
and k ~ ^m e /m^ ~ 0.07. Couplings to muons are mildly suppressed relative to the ones to taus, 
of the order \frn^Jm T ~ 0.24, which is not far from our understanding of the symbol (within 
order one). 

In Section |IV| we will encounter further viable patterns for A which are not limiting cases of 
either single lepton flavor Eq. ([Tj) or hierarchy Eq. (|8]). 


III. LFV PREDICTIONS IN TERMS OF R K 


The leptoquark framework with the A^ ansatz from Eq. ([8]) is very predictive and allows for 
correlations between rare 6-decays and LFV. We obtain 

B(B -» Kp ± e^ z ) ~ 3 • 10~ 8 k 2 


and 


~ 3 • 10 8 k 2 I -—— 

0.23 

~ 2 • 10 -8 k 2 I - —— 

0.23 


B(B -> ive ± r =F ) ~ 2 • 10~ 8 « 2 


B(B -> Kp ± r^) ~ 2 • 10 


-8 ( 1 - Rk \ ' 


0.23 


V 


B(p —>• eq) ~ 2 • 10 “ — 


_i2 k 2 (1 Rk 


B(t — > eq) ~ 4 • 10 — 


p 2 V 0.23 
-u k2 f^-Rx 


p 2 V 0.23 
-i4 1 ( 1 — Rk 


B(r M7 ) c 3 • 10- Q 23 

B(t —t prj) — 4 • 10“ n p 2 ( 1 - 


(13) 

(14) 

(15) 

(16) 

(17) 

(18) 
(19) 


For the amplitudes of the purely leptonic decays note the m T /m^ chiral enhancement of A(B & 
r + (//,e) _ ) over A(B S —> T~(p,e) + ) and A(B S -» p + p~). More general, neglecting phase space 

B(B S ->£+£'-) m 2 


B(B S -> £-£'+) mj 


( 20 ) 


We thank Franz Muheim for drawing our attention to future p —> e conversion experiments. 
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This relation follows from the left-handed lepton only hypothesis, a feature of the beyond SM 
(BSM) models considered here, see Appendix [A] for details. Eq. (201 can be used to test the lepton 
chirality in LFV processes. 

Furthermore we obtain 


B(B S —> p + e ) 
B(B S —> p + p~)sM 
B(B S ->■ r+e~) 
B(B S —> p + p~)sM 
B(B S ->• r+/i-) 
B(B S —»• p + p ~)sm 


~ 0.01 k 2 

~ 4 k 2 • ( 


/ 1 - B,k 
' V 0.23 
1 - R k \ 2 
0.23 ) 


— 4 • 


/ 1 — Rr 

V 0.23 


2 


2 


( 21 ) 

( 22 ) 

(23) 


where B(B S —> p + p ~)sm = (3.65 ±0.23) • 10 -9 [35]. The current bound B(B S —> p ± e T ) < 1.1 • 1CT 8 
at 90 % CL by LHCb [36] is hence at least two to three orders of magnitude away. 

In addition, for l ^ t', 

B(B -> Trgg?) B(B S ->■ Kt±e*) B(B^£±^) ~ ( Pd \ 2 ^ 

B{B -> A^±f±) ’ ^(5 AT#f±) ’ H(5 S -» ^±f±) V P ) ~ 1 ' 1 


Precise predictions for lepton flavor conserving decays we leave for future work. We stress that the 
calculation of the LFV hadron decays is much less complicated as contributions of quark loops or 
(gg)-resonances coupling to the electromagnetic current are absent. 

There are two extreme scenarios assuming Ao ~ 1: 

A) p = 0(1), corresponding to high mass leptoquark in the fewxlO TeV range, i.e., out of LHC 
reach. In this case k can be of order one, too, implying that B —> Kp ± e T and p —> eq could be 
just around the corner, cf. Table [I] Radiative T-decays are far away while B(B — > K(e, /Lt) = * = r =l= ) 
and B(t —> p,rj) are at least three orders of magnitude away from their respective current limits. 
B{B S -> t+h~) ~ KT 8 , B{B S -> T±e±) < 10“ 8 and B{B S -> p ± e*) < 10“ 10 . 

B) p -C 1, corresponding to light leptoquarks up to few TeV in mass and k< 1. As in scenario 

A), p —> eq can be just around the corner, B(B S —> t + p~) 10 8 and B(B —> Kp^r^) is three 

orders of magnitude away from its current limit. B —> K(t, p^e^, B s —> (p, and t —>• pp 

are strongly suppressed, but r —> pj can be found with the next round of experiments while r —> ey 
remains far away. 

C) Ao <C 1, corresponding to light leptoquarks up to few TeV in mass. This scenario arises for 
instance in the model in Appendix [B] where the leptoquark carries FN charge to increase powers 
of spurion insertions. The LFV pattern follows from the values of p , p d , k as dictated by the flavor 
model. 
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Concerning collider searches within pattern Eq. ([8j the leptoquark decays take place into second 
or third generation leptons and into third generation quarks. In scenario A), all combinations 
of quarks and leptons can arise at similar level except for those involving d-quarks, which are 
necessarily suppressed. 

In case of at least two more measurements rather than upper bounds on LFV in addition to 
Rk the parameters p and k can be determined, pinning down the flavor pattern of the leptoquark 
coupling matrix A further. It is conceivable that the latter is linked to the mechanism generating 
flavor for SM fermions, e.g., with flavor symmetries, which can be probed with the rare decay 
data. In Section |IV| we give examples of realistic flavor symmetries that give the single lepton 
flavor patterns Eq. ([7]), and special cases of the hierarchy pattern Eq. ([8]), but also further testable 
pattern. Predictions of leptoquark coupling patterns not covered by the ansatz Eq. ([8]) are given in 
Section IIVDI 


IV. FLAVOR SYMMETRIES 

In this section we illustrate how flavor symmetries generically control the shape of the leptoquark 
Yukawa coupling matrix A, and how there is in general a relation to Higgs Yukawa couplings 
controlled by the same symmetry. The reason for this is simply that both the leptoquark A 
and Higgs doublet (or doublets) couple to same fermions, whose generations transform as specific 
representations of the flavor symmetry. 

We use specific models to demonstrate this, and show how to obtain some special cases for A, 
including A^, A^ from Eq. ([t]) , and A^ from Eq. ©, but also new structrures. 

In terms of the quark index in A, the structure between rows depends on what kind of symmetry 
transformations we assign to either the Qi generations (for the LL leptoquark) or to the RH down- 
type quarks d? (for the RL leptoquark). In frameworks where the flavor symmetry distinguishes 
Q and d c , the LL and RL leptoquark models could be very different. Here, we focus mostly on 
frameworks of flavor symmetries explaining leptonic mixing, where this is not the case. If we assume 
all quarks are trivial singlets of whatever non-Abelian symmetry is included, one can still assign 
a FN charge to quarks (under U(1)fn) leading to Eq. © or similar. The FN mechanism relies 
on having lighter generations with larger charges in order to explain the mass hierarchies through 
added insertions of a 9 field. We choose without loss of generality to normalize the U(1)fn charge of 
6 to —1. If we assume that the leptoquark A is neutral under the U(1)fn, the leptoquark Yukawa 
couplings to lighter generations of quarks is simply suppressed by as many additional insertions 
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of 6 as the ones that appear for the respective masses, resulting in Eq. ©• We relegate further 
considerations of the U(1)fn charges of A to AppendixjBj in order to focus on the structure of the 
lepton couplings, where we consider non-Abelian flavor symmetries. 

Before we delve exclusively into lepton flavor models, it is interesting to consider a SUSY frame¬ 
work with an underlying 50(10) unified gauge group, where quark and leptons are linked. This is 
particularly restrictive in terms of the allowed flavor structures. At this level we will consider the 


flavor symmetry to be continuous ( SU(3)f , Section IVA). 

We then investigate in some detail a SUSY framework with discrete flavor symmetry A 4 in 
Sections |IVB| and |IV C| We focus on A 4 because it is a convenient framework to obtain the observed 
pattern of leptonic mixing, but also due to its relative simplicity as the smallest group with triplet 
representations. Appendix [C] contains a brief primer on A 4 . 

The reason for having SUSY in all the frameworks discussed in this section is twofold: first, 
SUSY keeps the gauge hierarchy problem under control, which is particularly relevant when going 
beyond the SM (with 50(10) or flavor symmetries broken at a high scale). Second, by holomorphy 
it allows one to separately align the different flavor symmetry breaking directions required in the 
respective models. Because of SUSY, we necessarily use two SU( 2) doublets, h u and h c j. 

All the frameworks also have a (7(1 )r R-symmetry, which is spontaneously broken to its Z 2 R 
subgroup acting like the MSSM R-parity. Another point that all the frameworks we consider have 
in common is the presence of an auxiliary Abelian flavor symmetry (either U(1)f, Z3 or Z4). As a 
point of notation, in each section the charge of a superfield cf) under the relevant auxiliary symmetry 
is denoted as {</>}. 


A. SU(3)f x U(1)f x U (1)r unified framework 

We discuss now a framework exemplified by the model presented in [ 8 j. Similar models were 
considered in [37| 138], showing how to obtain large reactor angle (# 13 ). One of the main features in 
the framework are that the 3 generations of each fermion flavor transform as triplets of SU( 3 ). 5 

In order to build flavor symmetry invariants, we add three flavons (each with 3 generations) 
named cf >|, <^ 3 and (the subscript numbers are labels, the upper indices are generation indices 
i = 1,2,3, anti-triplet indices of 5(7(3)). Through the details of the superpotential of the model, 
the flavons acquire vacuum expectation values (VEVs), breaking the flavor symmetry in specific 

5 Due to the underlying 50(10) unification, if the 3 generations of lepton 5(7(2) doublets transform as triplets, so 
must the generations of other SM fermions, as each generation is unified into a 16-plet of 50(10). 
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directions: 


(fa) = (0,0, a), < 023 } = (0,6, -b ), (0i 23 ) = (c, c, c) 


(25) 


with 5 ~ 0.20a, c ~ 0.206. This VEV hierarchy is required to be of order of the Cabibbo angle, but 
is also related to the hierarchy between the solar and atmospheric neutrinos |8j and even to the 
magnitude of the reactor angle 6 13 [38]. 

After 50(10) is broken to a left-right symmetric GUT containing 567(2)/, x 56/(2)# (before 
breaking to the SM), the left-handed fermions are referred as 0j and the conjugates of the right- 
handed fermions as 07 (where i,j = 1,2,3 are generation indices, triplet indices of 56/(3)), where 
-0 contains the 3 generations of Q, L and 0 C contains the 3 generations of u c , d c , e c and v c . The 
superfields containing the electroweak doublets we denote here as h u , hd- 

The Yukawa couplings appear from non-renormalizable terms which are generally of the type 
( 030 j)( 030 /)/i Uj rf and are controlled also by an auxiliary Abelian symmetry 6/(1)#. The VEV 


directions in Eq. (25) are responsible for giving hierarchical structures for all the charged fermions. 6 


We focus now on the leptoquark Yukawa structures A that occur naturally when leptoquarks 
are added to this framework. In order to study this, we need to keep in mind the transformation 
properties of the SM fermions, which are all triplets under 56/(3)#, and the charges of the flavons 
under 6/(1)#, which we denote in terms of curly brackets e.g. {0 3 } = 2, {02 3 } = 1 and {0i2 3 } = 3, 
while {h Uj d} = —4 and the SM fermions are neutral, {0} = {0 C } = 0. To predict A we simply need 
to select how A transforms under 56/(3)# x 6/(1)#. Here we discuss the simplest case where the 
leptoquark is an 56/(3)# singlet. In Appendix [D] we consider leptoquark triplets under 56/(3)#. 

Note also that, due to the underlying unification d c and Q both transform equally under the 
flavor symmetry so there is no difference in the A structures corresponding to a LL or RL leptoquark, 
and also that structures involving components from 0 ( Q , L ) and 0 C ( d c ) are symmetric due to 
their origin from the same 50(10) multiplet. If we were considering a unified gauge group that is 
not left-right symmetric then the Q and d c would generally transform differently under the flavor 
symmetry and the structures for the RL and LL leptoquark would be different. A prime example 
of this is 567(5), where lepton doublets L belong to the same GUT multiplet as d c , but Q is in a 
different GUT multiplet together with u c and e c . See Appendix [E| for some possibilities. 

If A is an SU(3)f singlet (as h u ^d), the flavor symmetry invariants shaping A are A Lj), 
which are very similar to those shaping the fermion masses. Which particular flavons couple to A 


6 Refer to [S] ES3 for the details related with getting viable mass ratios for the 2nd generation of fermions, such as 
m M /m s 0 1 at the GUT scale. 
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is determined by its U(1)f charge {A}: 


{A} = —2 leads to A((f> 2 3 Qi)((j) 2 ^Lj) which is a special limit of Eq. (J8J) where p = 1 and 
k = pd = 0, given in Eq. 


{A} = -3 leads to A + (4>hQi)(4i L j))- 

{A} = —4 leads to two contributions: the dominant one xA (<f>\Qi)(<j^Lj) and the subleading 


A (^23 < 2i)( < £l23- c 'j) + (^123 ( 3*)(^23- L l)j> shown together in Eq. <\27\. 
{A} = -5 leads to A ((^ 3 Qj)(^i23 L i) + • 


{A} = —6 leads to A^^sQiX^AA?)’ the democratic structure given in Eq. (28). 


As we intend to account for Rk, we show only the structures with non-zero entries simultaneously 
in X se , A be or simultaneously in A s/i , A^. 

The A flavor structure for {A} = — 2 reads 


At" 2 ! = A 0 b 2 


( 


0 0 0 

0 1 -1 

0 -1 1 


(26) 


which is a highly predictive limit of A^ Eq. ([8]) with the same exact magnitude on the 4 non-zero 
entries. As there is a single coupling we can absorb the order one coefficient by redefining Ao- 
Choosing {A} = —4 gives 


A'" 4 ] 


/ 


Aq be 


0 

1 


1 

2 


-1 \ 
0 


- 10^-2 


(27) 


which due to a > 6 > c has A couple predominantly to br. We kept x explicit to separate 
leading and non-leading contributions and redefine Ao to absorb the order one coupling of the 
subleading contribution, which, although suppressed by ~ bc/(xa 2 ) could in principle account for 
Rk- However, the simultaneous presence of the kaon constraint Eq. ( |11[ ) rules this out. 

The democratic pattern from {A} = —6 is 


At" 6 ! 


h i A 


A 0 c 2 


111 , 


(28) 


V 1 1 v 










13 


and very symmetric like A^ 2 !, but it preserves lepton universality so it can not account for Rr 7 ^ 1 . 

The only viable texture is therefore A^ 2 !. In fact the only obtainable textures if there are only 
these 3 flavon VEVs are the textures explored so far. {A} > 0 has no allowed couplings due to 
holomorphy, and for {A} < —6 the only possibilities repeat the existing patterns as the matrix 
structure is driven only be the flavons contracting with Qi and Lj. Even at higher order it is not 
possible with this held content to obtain linear combinations of the textures. The reason for this is 
clearer when considering a specific example, so we skip {A} = —7 as it has no allowed couplings, 
and take {A} = —8 where: 


^(^23^9)(^23^/1) ( e i?fc^3^23^123) ) 


(29) 


is invariant, with the same texture as A^ 2 ! but with overall U(l)p charge +6 added due to 
SU(3)f singlet 123 )- Even though U( l)_p would allow to swap flavon contractions to 

A((f>2 3 Qg)(4>^Lh)(eijk(j)2^(t > 23 ( l ) i23) aR d s ° on, a H except Eq. (29) automatically vanish due to the 
Levi-Civita tensor. 


B. A 4 x Z 3 x U(1)fn x U{ 1)r framework 

The first A4 framework we consider is of SUSY A4 x Z3 X U( 1 )fn models, initially proposed 
in [ 9 ] and with renormalizable UV completions m and in particular | 12 | which obtains non-zero 
reactor angle 6 13 in full agreement with current neutrino oscillation data. The FN mechanism m 
is implemented separately through U( 1 )fn, generating the hierarchy in the charged lepton masses 
without requiring small Yukawa couplings; it can easily be used to justify the quark hierarchies as 
described in Appendix [B] SUSY’s holomorphy together with the Z 3 separate the charged lepton 
sector and the neutrino sector. The A4 triplet flavons and <$> v acquire vacuum expectation values 
(VEVs) in special directions (1,0,0) and (1,1,1), respectively, shaping the leptonic mixing. The 
fields charged as 2 under the R-symmetry 1/(1 )r are the alignment fields and are responsible for 
giving the flavons their VEV directions. We do not discuss the details of SUSY breaking, but when 
it occurs U(l)u is broken to leave only a Z 2 R subgroup which distinguishes the SM fermions, i.e. 
it acts as R-parity. We refer to Table ITT] for details about the charge assignments. 

We have not listed the A4 representation of A in Table [TT] as it requires a more detailed discussion. 
As we have done for SU ( 3 )f, we focus on the case where the leptoquark is an A 4 singlet. We consider 
A 4 triplet leptoquarks in Appendix |D| 

For completeness, we include here very briefly the charged lepton Yukawa couplings in the 
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A 

Q 

u c 

d c 

v c 

L 

e c 


T C 

h d 

hu 

6 

< in 



€ 

C° 

a 4 

See text 

1 

1 

1 

3 

3 

1 

1" 

T 

1 

1 

1 

3 

3 

1 

T 

3 

3 

1 

^3 

2 

0 

0 

0 

2 

1 

2 

2 

2 

0 

0 

0 

0 

2 

2 

2 

0 

2 

2 

U(1 )fn 

0 

See App. 

B 


0 

0 

2 

1 

0 

0 

0 

-1 

0 

0 

0 

0 

0 

0 

0 

U( 1) R 

0 

1 

1 

1 

1 

1 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

2 

2 

2 

u( l)y 

1/3, 1/6 

i/6 

-2/3 1/3 

0 

-1/2 

+1 +1 +1 

-1/2 +1/2 

0 

0 

0 

0 

0 

0 

0 

0 


TABLE II: Field assignment within A 4 x Z3 yielding the single lepton patterns Eq. ([ 7 ]). Leptoquark as¬ 
signments are described in Section IV B We refer to Eq. |9| for quark FN charges, discussed in Appendix 


superpotential 

WL = ^ lL ^ eChd + ^° hd + x ^ lLrChd ’ ( 3 °) 

where A is a scale associated with the breaking of the flavor symmetry. Coupling to the (1,0,0) 
VEV results in a diagonal mass matrix for the charged leptons where we can identify L\ with e, L 2 
with /r and L 3 with r. 

For A, the hypercharge differs for the LL and RL leptoquark. The Yukawa coupling to lepto- 
quarks is associated with a renormalizable superpotential term that is either [d c ;ALj\ for SU(2) 
doublet A corresponding to Eq. (|2| - or, with a SU( 2) triplet A, one would have A^ [QiALj\ as 
a SUSY version of Eq. The quark generation index is i = d,s,b and the lepton generation 
index is j = e,n,T. We assume that all three generations of SU(2)l doublet and SU(2)l singlet 
quarks are trivial singlets of A 4 . We consider SU(2)l singlet down-type quarks as an A 4 triplet in 
Appendix [E] as this case arises in A 4 unified models [39l HOj with SU( 5). 

The columns of A are constrained because L is an A 4 triplet. If the leptoquark is an A 4 singlet, 
the renormalizable leptoquark superpotential terms are no longer A 4 invariant. A contribution to 
A appears at leading order (LO) from a non-renormalizable term where L contracts with an A 4 
triplet flavon. There are then three options for A to transform under A 4 , and 4 non-equivalent 
ways to build an A 4 invariant: 

1. A ~ 1. Then A isolates e, as M)Li] oc L\. 

2. A ~ 1". Then A isolates as oc L 2 . 

3. A ~ 1'. Then [0)Lj ] 77 A isolates r, as i]” OC L 3 . 

4. A couples equally to all lepton generations due to (</v)- 
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By coupling to the VEV (1,0,0), the different options for the leptoquark representation under A 4 
lead to a single non-vanishing column for A. Options 1. and 2. correspond respectively to explicit 
realization of the patterns A^ 6 ’^ in Eq. whereas option 3. gives SM-like Rk and is disfavored 
by the current LHCb measurement. Option 4., which occurs regardless of A ~ 1, 1", l', preserves 
lepton universality at LO but as discussed later has a myriad of LNU couplings already at NLO. 

For now we postpone option 4. and take the Z3 charge of A to be {A} = 2 (following notation 
used in other sections for Abelian charge). This choice means that at LO A couples to SM fermions 
only through 0/. 

This is a non-trivial result that we emphasize: using the same non-Abelian flavor symmetry A 4 
and VEVs that jointly predict viable lepton mixing with large angles, one can obtain automatically 
leptoquark flavor structures where LNU exists due to the isolation of a single lepton generation. 
This is consistent with the LNU in B to K decays, as shown in (3|, where the isolation of the e or 
/r generation was merely assumed. 

This isolation of lepton generation takes place at LO in a generic expansion parameter of (</>)/A, 
where (j) represents the flavons and A the scale of new physics associated with the breaking of A 4 . 
This can be identified more precisely in specific UV completions mm, where it may turn out 
the isolation is actually exact (due to e.g. missing messengers for NLO diagrams). Generically one 
may still associate next-to-leading order (NLO) effects to the presence of non-zero reactor angle 
013, which fixes the expansion parameter (cj ))/A ~ 0.2. 

With a sizable expansion parameter it is important to consider contributions beyond LO, i.e. 
terms featuring additional flavon insertions (of the triplet flavons or of a V of A 4 ). For the 
{A} = 2 case that produces lepton isolation, multiple insertions of 0/ involve the combination 
[cf>l(pijss which gives contributions to A that can be reabsorbed into the LO one, due to its effective 
(1, 0,0) VEV; on the other hand both <p u and carry non-trivial Z3 charge so they can only appear 
as multiples of 3 so the earliest contribution appears at NNNLO which is of the order ({(/>}/ A) 3 , 
hence sub-percent. We conclude that in this case, contributions beyond LO can only change the 
structure of A negligibly, and lepton isolation holds to a very good approximation. 

When considering effects beyond LO it is relevant to reconsider option 4., where {A} = 0 leads 
to LO lepton universality due to the coupling with <j> v . The Z3 always allows adding the neutral 
cj)l to this scenario, so there are NLO contributions involving [ 4 >i 4 >u] 3 a, 3 s- The effective VEVs of the 
combinations are respectively {4>\{ 0, <^,—<^)), {(j>\{2 (/>},,—4%,—($))■ For {A} = 0, A and 

[4>iL\ £'A are also allowed by Z3. LNU indeed arises at NLO, with 4 terms each with distinct LNU 
structures. Closer inspection of constraints in particular Rk Eq. ([5]) shows that deviations from 
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lepton universality at order ~ 0.2 in the entries A[ pK l 

( 


A 


NLO 


Ao 


Pd Pd Pd 
P P P 
1 1 1 


(31) 


are phenomenologically viable. We conclude that LNU at NLO gives a viable pattern for A. 

One could consider assigning quarks non-trivially under A 4 . however, in the present framework 
this leads to issues in obtaining viable quark masses and mixing. Instead we will explore this option 
further in a different framework, A 4 x Z 4 , in Section IIV C 1 


C. A 4 x Z 4 x U(1 )r framework 


The A 4 X Z 4 SUSY framework [TO] is another interesting A 4 framework. For its renormalizable 
UV completions see also P3, and [12] for considerations regarding large @13. 

We start with a brief comparison with the A4 x Z 3 framework discussed in the previous section: 
the first difference is that the VEV of the flavon triplet coupling to the charged leptons is now 
(4>l) oc (0,1,0). The second one is with respect with the FN mechanism and fields - neither U(1)fn 
nor 6 are present, with the charged lepton hierarchy being due to a field 6 '. 1 The last difference is 
that the sector neutral under Z 4 is the neutrino sector whereas with Z 3 , the charged lepton sector 
was neutral. This is particularly relevant due to e.g. the £ 7/ flavon neutral under Z 4 . 8 



A 

Q 

u c 

d c 

v c 

L 

e c 

P c 

T C 

hd 

hu 

& 

</U 

c 
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# 



a 4 

See text 

1 

1 

1 

3 

3 

1 

1 

1 

1 

1 

3 

3 

1 

T 

1 " 

3 

3 

1 

Z4 

2 

0 

3 

0 

2 

1 

0 

1 

2 

0 

1 

1 

0 

0 

1 

0 

2 
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0 


0 

1 

1 

1 

1 

1 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

2 

2 

2 

u( l)y 

1/3, 1/6 

1/6 

-2/3 

1/3 

0 

-1/2 

+1 

+1 

+1 

-1/2 

+1/2 

0 

0 

0 

0 

0 

0 

0 

0 


TABLE III: Field assignment within A 4 x Z 4 discussed in Section IV C It leads to patterns Eq. l8| with 


normal or inverted hierarchies regarding the lepton couplings. 


When adding a leptoquark to this framework one can again obtain LO A structures like Eq. 0 
by having the Z 4 charge of A be {A} = 2. At LO we have: 

1. A ~ 1. Then A isolates r as [WA oc L 3 . 


7 

8 


To avoid confusion with the of Section IV Bl we renamed this field to 9'. In I10H12I the same held is named £'. 
In this A 4 x Z 4 framework, the presence of allows viable #13 for a region of parameter space that is not fine-tuned, 
as discussed in detail in 1121. 





















2. A ~ 1". Then A isolates e as [((^}Tj] 7 oc L\. 

3. A ~ l'. Then [^L*] 77 A isolates p as [(<^)Lj] 77 oc L 2 . 
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They are similar to options l.,2.,3. of Section IVB modified slightly to account for the (0,1,0) 
direction of the VEV. 

More significant differences arise beyond LO, because the neutrino sector is neutral. Given that 
= 0, one is allowed to add it to the LO contribution thus constructing NLO contributions 


involving [(j>i(j>ii]3a,3s- This is similar to what we have seen in Section IVB although the effective 
VEVs of the combinations are now (<^ 2 (—0, <^)), (</i 2 (—</>*, 2<f> 2 , — <^)). 

In order to keep the model predictive, we assume now that the (j> v contributions to A are forbidden 
by a partial UV completion where beyond LO terms are only allowed through singlet insertions 
(this is a natural consequence if A 4 triplet messengers are absent [12]). In this scenario the only 
relevant flavon beyond LO is £ ,; , which is also neutral under the Z 4 . NLO contributions come from 
one insertion of and NNLO contributions from two insertions of The options for A are: 


1. A ~ 1 has LO coupling to r, NLO coupling to electron, [<^Lj] 77 £ // £ ,/ A NNLO 

coupling to p. 

2. A ~ l" has LO coupling to electron, \4>]Li ] 77 £ 7/ A NLO coupling to p. [0jLj] £"£,"A NNLO 
coupling to r. 


3. A ~ 1' has LO coupling to p, [^Lj] A NLO coupling to r, [^Lj] 7 A NNLO coupling 
to electron. 


We illustrate these in matrix form, defining k" = (£ /; )/A ~ 0.2: 



^ n"p d K tr2 p d p d ^ 


^ Pd K"Pd K" 2 Pd ^ 


^ n" 2 p d Pd n"pd ^ 

aM ~ Aq 

k" P k" 2 P p 

, A^ ~ Aq 

P K "p k" 2 P 

, A^ ~ Aq 

k" 2 p p n"p 


v k" k" 2 1 j 


i 1 K" K" 2 j 


v k" 2 1 re" , 


(32) 

where we also parametrized the quark suppression factors as in Eq. Q for a more direct comparison. 
The first pattern, AW cannot simultaneously accommodate Rk and i?, s -mixing constraints unless 
M < few TeV. The second one, corresponding to an inverted hierarchy with leptons flavor ordering 


of leptoquark couplings in opposite way as the one to the Higgs, and the third one, normal hierarchy 


(for electrons) are both viable. Phenomenological implications are summarized in Section IV D 
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To conclude, within A 4 X Z 4 it is possible to obtain special versions of Eq. 0. Note also that, 
if instead of we had considered a Z 4 neutral V flavon - which would be difficult to distinguish in 
terms of lepton mixing angles |12] - for the same LO contribution, replacing the 1” with the 1’ leads 
to swapping which lepton generation is coupling at NLO and at NNLO, i.e, effectively swapping in 


each matrix in Eq. (32) the k" with the k " 2 terms. The resulting A structures are different and can 
be tested experimentally. As a corollary of that, by having simultaneously both V and one can 
obtain LO in one lepton generation and at NLO the other two, e.g. LO coupling to /i, NLO to e 
and to r. The variant with LO coupling to t is viable only for leptoquark masses not exceeding a 
few TeV due to Eq. 0 and Eq. 0. 


1. Quarks non-trivial under A 4 

An interesting option that can be considered is to make different generations of quarks transform 
as different non-trivial singlets of A 4 . Consider option 1., 2., 3. of the A 4 x Z 4 framework. If all 
quark generations are the same (non-trivial) A 4 singlet, it just shifts which lepton generation is 
isolated at LO by each leptoquark choice. For instance, if Qi are all 1/, then A ~ 1" is the 
leptoquark that would now couple to r, whereas A ~ 1 which previously isolated r would now 
couple instead to [i. 

On the other hand new flavor pattern arise when different generations of quarks are transforming 
under different A 4 singlets. We note that some care is necessary as this possibility may lead to 
unwanted implications for the Yukawa couplings with the Higgs and a CKM matrix that is not 
viable. Indeed this is generally the case for the A 4 X Z 3 framework. 

What happens with different generations of quarks assigned as different A 4 singlets is that, 
depending on the quark generation, the leptoquark couples to different columns at each order. As 
an illustration of this, and neglecting the quark FN charges for simplicity, take A ~ V , d\ 3 still as 
trivial singlets but d\ ~ . We have at LO: 

A (.X2<i2 + x sd 3 ) + A [x\d\) (33) 

corresponding to: 

f 0 0 x\ \ 


A LO = A 0 


0 12 0 
yO x 3 0 j 


(34) 





19 


Such patterns generalize lepton flavor isolation Eq. ([7]), and at the same time are no longer special 
limits of Eq. (|8]) . The pattern Eq. (34) could be distinguished from the others for instance through 
the tree-level decays of the leptoquark in colliders. 

Viable masses for quarks can be obtained for this type of d c assignments under A 4 as long as 
there are flavons allowing invariants populating the respective columns of the Yukawa coupling 
to the Higgs, h^QiY^d^. This is not the case for the A 4 X Z 3 framework, but in the A 4 X Z 4 
framework considered in this section the Z 4 neutral flavon allows it. For the example considered 
above where d\ ~ l ", terms populating Yf- with j = 1 can be obtained through a double insertion 
of the flavon £" (note that if there is a FN charge of d\ it should be amended with the additional 
flavon insertions in mind cf. Appendix 0- For this reason, choices leading to Q 3 u 3 non-trivial 
under A 4 are unappealing from a theoretical standpoint, as having the top quark mass arise from 
a renormalizable term is preferable. 

In general, constructions of this type yield patterns for \ LO where a single quark flavor can 
couple to a single lepton flavor. Due to Rk, either electrons or muons have to couple to both s and 


b. Therefore, the following structures in addition to Eq. (34) are phenomenologically viable: 

( xl 0 0 \ / 0 0 si \ / 0 X! 0 

12 0 0 j 

^13 0 0 J 


\ 

0 12 0 
\ 0 x 3 0 J 


\ 

X 2 0 0 

\ X 3 0 0 J 


(35) 


All of these effective two-flavor patterns predict LFV involving two lepton flavors only at LO, but 
receive respective NLO (and NNLO) contributions which are similarly shifted in the first row and 
are also themselves a two-flavor pattern. 

However, in specific UV completions it is possible to forbid NLO and NNLO contributions to 
A, while allowing the desired LO contributions and the necessary Yukawa couplings. 9 In fact, 
allowing only a minimal set of messengers we can obtain A matrices like Eq. ( [34] ) with x\ = 0, 
corresponding to the truly minimal patterns from |3] that explain LNU data Eq. ([I]). 


D. Summary of flavor symmetry frameworks and predictions 


We summarize in Table IV the features of the SU(3)f and A 4 frameworks discussed in the 


previous sections, for a single leptoquark. 


9 A relatively simple possibility for RL leptoquarks is to forbid A 4 triplet and SU(2) singlet messengers. We have 
checked in this case the absence of pattern-changing beyond-LO diagrams, while the required diagrams for RL 
leptoquark Yukawa couplings and the hd Yukawa to quarks and charged leptons are mediated by SU(2) doublets. 
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A structure 


Eq. (81 with zeros, i.e. Eq. (261, 2 1 


Eq. (7), A^l electrons only 
Eq. (7), A^l muons only 


Eq. (8) with k~ 1, i.e. Eq. (311, A NLO 


Eq. (8), normal hierarchy, i.e. Eq. |32|A^ 1 
Eq. ( 81 , inverted hierarchy, i.e. Eq. ^32} A^ 1 1 


Eq. (34) , Eq. (351 two-flavor 


symmetry 

flavons 

A assignment 

SU{3 ) f x U(1) F 

<^23> = (0,6,-6) 

{A} = —2 

A4 X Z3 

II 

0 

1, {A} = 2 

A4 X Z3 

1? 

II 

'T 

0 

0 

1”, {A} = 2 

^ 4.4 x 

(<l>v) = (W,W,'W) 

l x , {A} = 0 

A4 X Z4 

<0,) = (O,«,O), C 

r, {A} = 2 

A4 X Z4 

II 

e 

«nrv 

1”, {A} = 2 

A4 X Z4 

= (0, u, 0) 

1’ , {A} = 2 


TABLE IV: Summary of viable A structures and how they can be obtained through a flavor symmetry. In 


all cases except Eq. (261 and Eq. (34) the quarks transform trivially under the flavor symmetry. 

The phenomenology of the patterns following from the hierarchy ansatz Eq. (|8|) has been detailed 


in Section III The two-lepton flavor patterns Eq. (34) and Eq. (35) discussed in Section IV C 1 that 


arose newly in the analysis of the flavor symmetries, have predictions similar to the single lepton 
flavor pattern Eq. Q: they successfully explain Rk by either BSM in b —> se + e~ or b —> s/r + ^ - 
transitions, however LFV signals can appear in addition in either e/r, er or fir related to b —> d 


or s —> d transitions. For the example in Eq. (34): B —> 7 tt/x, B —» r/r and r —> K^/i, but 


none involving electrons. The sub-case for x\ = 0, that is no couplings to d-quarks, can also be 
obtained, corresponding to the minimal benchmark patterns presented in [3j, to which we refer 
for their phenomenology. They have no LFV. Implications include effects in b —> (d,s)vv decays, 
relevant for the Belle II experiment. Correlations between B —> Kvu and B —)• K*vv have been 
worked out in EJ 

The other genuine non-Abelian type of A structure is the inverted hierarchy Eq. (32), A^ 1 1, in 


which the leptoquarks couplings do not follow the ones of the leptons to the Higgs. This pattern 
predicts LFV in processes involving ep, er and /jt, which can be read off from the estimates given 


in Section III using k = n", n = n" 2 and k = n //3 , respectively. The normal hierarchy pattern, A^ 1 1 
of Eq. (32), predicts LFV in processes involving e/j ,, er and fir corresponding to k = k" 2 , k = k" 3 


and k = k", respectively. In these cases, k" = 0(di3). 

While it appears that it is possible to obtain in each flavor model pattern for RL and LL 
leptoquarks alike, we stress that this is a feature of the specific frameworks considered where Q 
and d c transform equally under the flavor symmetry, either due to underlying unification or for the 
sake of simplicity. 
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V. RARE HIGGS DECAYS 


In this section we discuss leptoquark effects in decays of the Higgs into fermion anti-fermion. 
Such decays have received recent interest with the advent of such a particle and its various decay 
modes being stringent test of the flavor sector of the SM and beyond, e.g., recently, 021 S3]. 

Leptoquark contributions to Higgs decays are induced at one-loop as exemplified in Figure [lj 
The amplitude is proportional to a renormalizable term At A h)h, whose coefficient in general is 
model-dependent. After electroweak symmetry breaking a coupling of the Higgs to two leptoquarks 
is induced at order vhAA*, where v ~ 174 GeV sets the electroweak scale. The corrections to 
flavor-diagonal modes h —> It are hence parametrically given as 


tyi ~ N y2 

yi c M 2 16-7T 2 ’ 


q = d,s,b, 


(36) 


(37) 


where N c = 3 denotes the number of colors. Analogously for h —> qq\ 

5y q v 2 \\ q i\ 2 B _ 

y q ~ M 2 167T 2 ’ e, ^’ T ' 

Contributions from further one-loop diagrams involving the Higgs-Yukawa couplings exist but are 
suppressed further by at least yt, T <S 1, and contributions with yt don’t bring in additional en¬ 
hancements either. The relative corrections 5y/y to the diagonal Higgs Yukawas therefore do not 
exceed the 1CP 2 level. 


A 

/ 


/ 


-< 


v 


\ 

\ 

A \l 


x* 


q 

Kjf 


£ 


m 


FIG. 1: 1-loop contribution to h —» ££' from A. 

Spurion analysis shows when LFV in the Higgs coupling liXeR is induced, cf. Eq. Q, 

X = AY e + A Y e + ... , (38) 

hence the leptoquark couplings, in general, induce LFV couplings. A single Higgs doublet as in the 
SM suffices. Note that single lepton-flavor patterns as in Eq. ([7]) fail to induce Xij for i ^ j. 
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Leptoquark contributions to the H-term in the spurion expansion arise at order 


yw 

ye' 


N r 


V 2 A*^A, 


q e A q(-' 


c M 2 16tt 2 ’ 

yin is suppressed relative to yw by yi/yi' for mi < my. 


q = d,s,b. 


(39) 


The coupling in Eq. (39) is the same which drives £' —> £7. We find that, using the limits from 


Table |TJ LFV Higgs effects are limited as 


V JN < i 0 -3 ! EL < iq-3 ^ < 10 

y tJ . 


yer 


lJe/i 


(40) 

2/r Vt 

Flavor models of course predict patterns for these couplings and relate h —> f!i decays. E.g., in 
terms of the parametrization Eq. (pi): 


Her Uefi 

- ~ K , —— ~ K . 


(41) 

V/J.T yfir 

If k is large, this implies that all LFV branching ratios are strongly suppressed due to the suppression 
of h —> ye. If k is small, the h —> r/r branching ratio is much larger than the ones involving electrons. 
Different patterns follow from other flavor structures of A summarized in Section |IVD 

The CMS data on the h —> ry branching fraction [|7j imply a sizable off-diagonal Yukawa- 
coupling to the SM-like Higgs ^/\y T/J ,\ 2 + ly^A 2 = (2.6 ± 0.6) • 10^ 3 , about a third of the coupling 
to the taus itself, y T — 0.01 tan/3. Such a sizable effect, taken at face value, exceeds the leptoquark 


estimates Eq. (40) and points to mechanisms beyond perturbative loops. An example is given by 

mm- 


VI. CONCLUSIONS 

The bottom-up BSM scenario with leptoquarks considered here provides further evidence of 
how new physics can help to learn about flavor |46j : if the anomalous signatures in the FCNCs 
Eq. ([!]) hold they offer unique possibilities to probe the mechanism of flavor, such as the type and 
charge assignment of a flavor symmetry whose imprints on SM matter were the only experimental 
information previously available. 

By adding a leptoquark to existing flavor symmetry frameworks we obtained sample scenarios 
where the flavor symmetries are simultaneously responsible for LNU ( Rk ) and LFV, drawing in 
each case relations to the lepton mixing angles and charged fermion hierarchies predicted in these 
models. 

The predictions for LFV lepton and ^-decays are copious, and we summarize them in Sections 


III and IV D The LFV branching ratios, whose size is driven by (1 — Rk) 2 are in part accessible to 
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near-term future experiments, including MEG, LHCb and Belle II, the details of which depend on 
the specific flavor pattern of the leptoquark coupling A, single lepton flavor Eq. ([7]), the hierarchy 
pattern Eq. pi), or further patterns which follow from the non-Abelian nature of the underlying 


flavor symmetry. Viable and realistic ones are summarized in Table IV They can be distinguished 
phenomenologically. Further data could not only be a discovery of BSM, but also provide clues about 
flavor. Already current data rules out many choices of flavor group assignments, as demonstrated 
in Section 113 

Indeed if certain patterns of leptoquark coupling are observed in the future, they constitute 
strong hints of the presence of non-Abelian flavor symmetries. From the examples studied in this 
paper we single out three patterns. One is predicted in the unified SU(3)f framework, and is highly 


symmetrical with four entries of the same magnitude, Eq. (26). Alternatively if lepton isolation, 
either in electrons on muons, Eq. ([7]), is observed to below percent level, that points to models such 
as our A4 x Z3 which only allow deviations from isolation at the NNNLO - or, given an appropriate 
UV completion, to models like our A 4 X Z 4 . Mimicking these patterns without using non-Abelian 
symmetries would require rather unnatural leptoquark couplings. 

If LNU and LFV in rare 6-decays is observed, it allows to completely identify the chirality- 
structure of the leptoquark-lepton-quark couplings: comparison of Rk with related non-universality 
tests into others strange final states, such as K*,X S ,... allows to probe for right-handed quark 
currents m Specifically, the models Eq. Q and Eq. ([ 3 ]) can be distinguished. Pinning down 
the chirality of the leptons is possible by comparing the LFV branching ratios B^ —> l + l'~ with 
B^s) —> which pick up different lepton mass factors. This method of diagnosing quark and 


lepton chirality, of course, applies model-independently for the low energy theory Eq. (Al), i.e., is 
not limited to leptoquark mediated FCNCs. 

The current hint for BSM in \AB\ = 1 transitions in Rk provides an anchor for fixing BSM 
scales when interpreted as a “signal” and combined with |AB| = 2 (IL,-mixing) bounds. The mass 
range of the leptoquarks is determined to be right above search limits around O(TeV) and below 
(D(50)TeV [3]. In pp collisions the leptoquarks will be pair-produced and decay in our frameworks 
predominantly into second or third generation leptons and into third generation quarks. 

Note added: During the completion of this work a paper (48), appeared, where related lepto¬ 
quark effects to h —> t(! were considered. 
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Appendix A: Leptoquarks in 6-decays 


This appendix explains how to obtain predictions for 6-decays from the leptoquark interactions 
Eq. Q and Eq. (|3j), corresponding to A(3, 2)x/g and A(3, 3 )_i/ 3, respectively. After integrating-out 
and fierzing —> s£ + £' transitions with £ = £' and £ ^ £' are induced. Employing the common 

|AR| = |A«Sj = 1 effective Hamiltonian 

n eS = -4 C -Jv tb v*2 E c *°' ’ ( A1 ) 

V i 

where a e , V t j and G F denote the fine structure constant, the CKM matrix elements and Fermi’s 
constant, respectively, the semileptonic operators 


09 = s llx P L b ff'fl , 010 = s^P L b £'^ 75 £, (A2) 

0' 9 = s^P R b£'r£, O' 10 = s llx P R b£'^ 5 £, (A3) 

receive leptoquark contributions as |,3] 

for A(3, 3)_ 1/3 , 

for A(3, 2) 1/6 . (A4) 

Note that + c [q = 0 shows that the leptons are left-handed. 

The amplitude for B s —> £ + £'~ decays is calculated as 




0 9 = -0io = 


0/ _ 
°10 — 


Vr \* e \u V2 
a e V tb V t * s 2M 2 G f 

_£»/ _ 71 ^se^be 1 

9 “ a e V th V* 4M 2 G f 


(0\s lfl P L/R b\B s ) ■ u(t 


0$V + 0 


l oV75 


v(£ + ) ex mt,{c [2 + 0 f) + m e (C§ - C^) 


(A5) 


where we used that the hadronic matrix element is proportional to the B s mesons’ four-momentum, 
which equals the sum of the leptons’ four-momenta q and then applied the equations of motion 
for particle u and anti-particle spinors v. For same flavor leptons therefore the vector current 


contribution vanishes and using Eq. (A4) the r.h.s. of Eq. (A5) reads 2rn^C^ 9 = —2mgCg\ with 








25 


proportionality given by the anti-particle mass. Right-handed leptons, corresponding to —c[|j = 
0 would give instead the mass of the negatively charged lepton as chiral factor. 


The model-independent framework Eq. (Al), Eq. (A4) allows to compute leptoquark effects in 
B -> pp decays in a straightforward way. Correlations of the global fits and Rk in such models 
have been discussed recently in mm- A detailed exploration including the very recent, preliminary 
3fb _1 data by LHCb on B —> K*pp angular observables [19] taking into account SM uncertainties 
is beyond the scope of our work. 


Appendix B: Quarks and Froggatt-Nielsen 


Relative suppression between rows of A is naturally obtained by assigning different U(1)fn 
charges to the generations of quarks. The outcome depends also on the charge of the leptoquark. 
In this appendix curly brackets denote the U(1)fn charge of a field. 

Explicitly, take the leptoquark SU( 2) triplet, {A} = 0 and {Qi, 2 , 3 } = 3,2,0 respectively 
The leptoquark couplings would be: 


(x d Qi Ae 3 + XsQoAe 2 + x b Q 3 A) i , 


(Bl) 


where e = (0)/A ~ 0.2. Xd, s ,b are order one dimensionless couplings. The outcome is an hierarchy 
between rows of the leptoquark couplings X d e < X s e < as e.g. \ d i/\ s t = yje, irrespective of 
lepton flavor. One can easily generalize this for alternative assignments of FN charges, such as 
{ 1 , 2 , 3 } = 4,2,0 |21j . and similarly for the leptoquark SJJ{ 2) doublet, where {d\ 23 } = 2,1,0 or 
{du 3 } = 1 , 1,0 m Values for suppression factors of strange quarks versus 6-quarks, p, and down 
quarks versus 6-quarks, p d , can be read-off and are compiled in Eq. 0. In this framework p d /p < 1, 


consistent with the bounds Eq. (12). 


Choosing non-zero values of FN-charges for the leptoquark significantly impacts A and the 
resulting phenomenology, which forces {A} to be within 0 and ~ 3: {A} > 0 adds an overall 
suppression to all leptoquark couplings of which, while keeping p d /p unchanged, leads to 

suppression of Ao- This requires lighter leptoquarks to explain Rk- This case is discussed in 
scenario C of Section III Too large charges {A} > 3 cannot accommodate simultaneously Eq. 0 
and direct search limits M > 1 TeV. Within SUSY negative charges are in conflict with Eq. 0, 
because due to holomorphy, {A} < 0 forbids couplings to 6 -quarks. To see this take {A} = — 1 and 
compare to {Q 3 } or {d§}, which are zero to account for the mass of the b and top quarks. 
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Appendix C: A few A4 details 


A4 has 4 distinct irreducible representations, three of them are 1 -dimensional, i.e. singlets, and 
one of them is 3 -dimensional, i.e. the A4 triplet. The trivial singlet we label 1 , and it transforms 
trivially under A4. The (non-trivial) singlets l 7 and 1 " are conjugate to one another and they 
transform under a specific A4 generator, T, by getting multiplied respectively by tc 2 and ui (oj = 
e « 27 r/ 3 , w ith u) 3 = 1 ). The product of A4 singlets has (1 x 1 '), (1 X 1 "), (l 7 x 1 ") transforming as 
l 7 , l 77 , and 1 respectively. The action of the group on triplets is represented by 3 X 3 matrices. 
Consider specific triplets A = (01,02,03), B = (61,62, ^3)- Under generator T, which is diagonal in 
the basis we are considering, we have TA= (ai, ca 2 02, coa 3 ) (same for B). We use the conventions 
in mm, and square brackets to indicate A4 products: 


[AB] =(ai&i + 02^3 + 03^2) ~ 1 , 
[AB]' =(0162 + a-2h + a 3 b 3 ) ~ l 7 , 
[AB]” ={a\b 3 + a 2 6 2 + a 3 b\) ~ l 77 . 


It is also possible to construct a symmetric (s) and an anti-symmetric (o) triplet: 


/ 


[AB]s = 3 


2ai6l - 02^3 — 03^2 


\ 


20363 - Oi&2 - 02^1 

\ 20262 - 03^1 - a \ b 3 


( 


[AB} a = ~ 


0263 - a 3 b -2 


\ 


Oi6 2 - 02^1 
\ fl 3 &l - 01^3 / 


(Cl) 

(C 2 ) 

(C 3 ) 


(C 4 ) 


Appendix D: (Anti-)triplet A 


Within the frameworks we considered in Section IV it also possible to consider 3 generations 
of leptoquarks transforming as a representation of the flavor symmetry. This is more intricate 
than the single leptoquark scenarios we focused on, as there are three A structures (one for each 
generation of the A multiplet), and the flavor structure of the mass matrix for the A generations 
is also relevant. In the frameworks considered here A only acquire masses after the R-symmetry is 
broken (the leptoquark mass terms behave similarly to the /i-temi /r/r. u /i^). As we will show, in some 
cases there are holomorphic A bi-linears that could couple to whatever superfield is responsible for 
breaking the R-symmetry. 
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1. SU(3)f 


In the SU( 3 )f framework one can have anti-triplet A* (i.e. 3 generations of leptoquark trans¬ 
forming like the 4 > l ) or instead a triplet A.; (i.e. 3 generations of leptoquark transforming like the 
SM fermions). Discussing full models is beyond the scope of the present paper, but we illustrate 
some structures that can arise. 

For triplet A,;, the invariants can be quite different from those discussed in section IV A If 
{A} = 0 there is an invariant not involving the flavons e* jfc A iQjL^. This is a purely anti-symmetric 
structure for each A*: 


/ 


A 1 = A 0 


0 0 0 

0 0 1 

0-10 


\ 


/ 


A 2 = A 0 


0 0-1 

0 0 0 

1 0 0 


\ 


/ 


A 3 = A 0 


0 1 0 

-10 0 
0 0 0 


\ 


(Dl) 


For this charge assignment there are no holomorphic leptoquark bi-linears to study, but we can 
nevertheless conclude that this type of structure can not account for Rk regardless of the mass 
eigenstates. 

If {A} ^ 0 , invariants can arise by contracting to 3 flavons, such as 


A 1 = A 0 


0 0 0 
0 0 0 
0 0 0 


A 2 = -A 3 = A o & 3 


0 0 0 

0 1 -1 

0 -1 1 


(D 2 ) 


which occurs for {A} = — 3 , where Ai decouples and A2,3 couple with A 2,3 ~ A^ 2 1 , Eq. ( 26 ). For 
this charge assignment, there is a holomoprhic bi-linear (</>|Aj)(</>2 3 A j) + (i^Aj^^gA j). It could 
give rise to a mass contribution after the R-symmetry is broken: 


Ma 


m/s.ab 


^0 0 0 
0 0 1 



V ° 1 - 2 / 


(D 3 ) 


Diagonalizing this mass matrix reveals Ai has vanishing mass prior to soft terms, so we expect it to 
remain as the lightest leptoquark - however it also has no coupling to SM fermions. The next lightest 
mass eigenstate with mass mass ( — 1 + y/2)m^ab, and the heaviest with mass (1 + y/2)m/±cib , are 
respectively dominantly A2 (with some A3) and the orthogonal combination. As both A2,3 couple 
to SM fermions with the A^ 2 ^ texture, this option reduces to a more involved version of the singlet 
leptoquark discussed before. 
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For anti-triplet A*, the invariants require only a single flavon contraction (contrasting to fermion 
mass structures that required 2 , and the anti-triplet which as we have seen requires either 0 or 3 ). 
An example would be for {A} = — 1 , where we have (A l Qi)((fx 23 Lj) + ((j>2 3 Qi)(A J Lj): 


Ai — Aq b 


0 1 -1 

1 0 0 

-10 0 


A 2 — Aq6 


0 0 0 
0 2-1 
0-10 


A3 — Aq b 


0 0 0 

0 0 1 

0 1 -2 


(D 4 ) 


There is no contribution from holomorphic bi-linears, as e^k^^Ai A k vanishes. We can still con¬ 
clude that these structures would allow for Rk 7^ 1 as long there are sufficiently light eigenstates 
containing A 1,2 . 


2. A 4 


If we have 3 generations of leptoquarks as an A4 triplet A ? ; there can be invariant contractions 
[LA] = L\A\ + L2A3 + L3A2, if under Z3, {A} = 2 ; within A4 X Z4 the choice would be {A} = 3 . 
Each generation has its own A* matrix. The LO-structures are: 


A 1 


Aq 


( x d 0 0^ 

i s 0 0 

\x b 0 0 ) 


A 2 


Ao 


^0 0 
0 0 x s 
\ 0 0 x b I 


A 3 


Ao 


' 0 x d 0^ 
0 x s 0 
V 0 x b 0 ) 


(D 5 ) 


The holomorphic bi-linears in A allowed by A4 x Z3 are x^ [AA], x'^' [AA] /7 and x u [(f> v AA\. This 
could lead to a mass structure 


M a 




(\ 0 o\ 

0 0 1 

V° 1 V 


+ x^t 


(0 0 l\ 

0 1 0 

V 10 V 


1 

+ 3 *. 


/ 2 -1 -l\ 
-1 2 -1 
-1 -1 2 


(D6) 


where we absorbed the magnitudes of the respective VEVs into the coefficients x. This is similar 
to the Majorana neutrino structure in this framework (see [ 12 ]). With eigenvalue + x^, ( 1 , 1 , 1 ) 
is an eigenstate of this structure which, when compared to the respective A*, would preserve lepton 
universality. In the limit (£') = 0 we have effectively x'^ = 0 and the other two eigenstates are 
also independent of the free parameters: they would be (2, —1, —1) and (0,1, —1), respectively with 
eigenstates x u + x% and x v — x^. These two leptoquark eigenstates could mediate LNU couplings. 
As it is (£') ^ 0 that generates non-vanishing reactor angle in this A4 x Z3 framework, in realistic 
regions of parameter space the (2, —1, —1) and (0,1, —1) directions are no longer eigenstates of the 
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leptoquarks. Rk 7^ 1 is still possible and an interesting situation arises where the perturbations 
away from the (2, —1, —1) and (0,1, —1) directions to the leptoquark mass eigenstates are directly 
related to $13 7I 0 and the required perturbations of the neutrino eigenstates, appearing in the 
leptonic mixing matrix (recall in this basis the charged leptons are diagonal). 

Studying NLO corrections for the A\ triplet A t . which requires correcting both the A* and the 
assumed structures for Ma, is beyond the scope of the present work. 


Appendix E: d c as A4 triplet 


We consider here a situation where Q remains as singlets under the flavor symmetry but d c is, 
like L, an A4 triplet. This situation arises naturally in A4 unified models of lepton mixing [ 39 l; 40 j 
with SU( 5 ). Viable down quark masses can be obtained with LO contractions to the 4 >i flavon, with 
NLO corrections from non-trivial singlet flavons and the up sector enabling viable CKM mixing. 

If d c ~ 3 one can have a RL leptoquark coupling with no flavon, [d c L\ A^A, which depending on 
A ~ 1 , V , 1 " gives respectively A^ structures: 


A [1 = Ac 


(\ 0 o\ 

0 0 1 

V° 1 V 


A [1 = A 0 


/o 0 l\ 

0 1 0 

v 10 °y 


AH = A 0 


(0 1 o\ 
1 0 0 

V ° 0 V 


(El) 


This structure would appear for {A} = 2 within A4 x Z 3 and for {A} = 3 within A4 X Z4. 

Another possibility would be to couple to 4 >i, [d c L\ a}S (/)i\^ A. For a (rt, 0 , 0 ) VEV and 
A ~ 1 , T, l", or alternatively for a ( 0 ,ii, 0 ) VEV and A ~ 1 // , 1 , 1 / , these structures correspond 
respectively to A^: 


^2 a s 

0 

0 ^ 


( 0 

T 

0 ^ 


/ 

0 

0 

da T 

0 

0 

&a H - 0's 

1 

& a H - 0'S 

0 

0 

5 


0 

2 CLs 

0 

V° 

& a “1“ 0'S 

0 J 


\ 0 

0 

2 CLs ) 


\a a 

+ CL S 

0 

/ 


One of these A^ can appear alone for the A4 x Z4 framework and {A} = 1 , as { 4 >i} = 2 . In the 
A4 x Z3 framework {<(>;} = 0 so the A^ structure appears for {A} = 2 , i.e. together with the 
respective A^ structure for each choice A ~ 1 , 1 ', 1 " - the choice V having A se ,be 7^ 0 and the choice 
1" having A7^ 0, these are two more candidates that can justify Rk 7^ 1. 

The final possibility is to couple to <f > v , with ( 1 , 1 , 1 ) VEV, through [d c L\ atS <f) v A^A. A ~ 1 , 1 ', 1 " 
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correspond to 


AM = 


0 +1 -1 
-1 0 +1 
1 -1 0 


( 2 -1 -l\ 


+ ds 


-1 2 -1 

-1 -1 2 


(E3) 


In the A 4 x Z 4 framework {4>i/} = 0 so the A^ structure appears for {A} = 3, i.e. together with 
the respective A^ structure for each choice A ~ 1, 1', 1". A^ can appear for the A 4 X Z3 framework 
and {A} = 0, as {<f> v } = 2, but this also allows A to couple to Effectively this means that in 
both frameworks one must combine A^ with one of the three A^. 
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